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' Abstract 

m 

^ ■ Proton decay is examined wilhin the fnmiowork of certain S()( KJi models 

P< which have low dimension representations compatible with the recently found 

constraints coming from string theory, and which use the Dimopoulos-Wilczek 

Si 



mechanism to solve the doublet-triplet splitting problem. It is found that the 
mass parameter that controls proton decay rate is intimately related to the low 
energy data via the coupling unification condition, and that a suppression of 
proton decay is achieved at the cost of large threshold effects. Furthermore, 
in cases where there are states with intermediate masses, the threshold ef- 
fects further enhance proton decay. The experimental bound then severely 
constrains the parameter space of these models. Some possible solutions are 
suggested. 
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SO(10) is an attractive gauge group for grand unification. It may be thought of as 
arising via gauging the global B-L symmetry of SU(5) models. In SO (10) models, all the 
quarks and leptons for each generation are unified into a single representation, 16, and since 
SO (10) observes left-right symmetry, anomaly cancellation is automatic. Also, the Yukawa 
couplings for the members of a family (e.g., top, bottom, and r for the third generation) 
must be the same at the GUT scale and this has allowed a successful prediction for the mass 
of the top quark, provided that the value of tan/3 is large, ~ 50. Another nice feature of 
SO(10) is the existence of a way to avoid the doublet-triplet splitting problem via the so 
called "Dimopoulos-Wilczek mechanism." Jl].^] 

A drawback to SO(10) models in the past has been the largeness of the particle content, 
due to the use of large representations such as 126 and 210. From the viewpoint that a 
grand unified theory (GUT) is an effective field theory of some Planck scale theory which 
is perturbative below the Planck scale (e.g., a string model), it is desirable not to have a 
large number of chiral multiplets since these then would drive the unified coupling constant 
upwards, exceeding the perturbative bound before reaching the Planck scale. This desire 
is consistent with recent results in free field heterotic string models 0, where it has been 
shown that no massless SO(10) representations larger than 54 can arise, regardless of the 
affine level at which SO(10) is realized. 

There are now several SO(10) models where the Dimopoulos-Wilczek mechanism is em- 
ployed while using only small particle content ||,|5| . Here, we will refer to these as "minimal" 
SO(10) models. In this letter, we present two new results that are relevant to proton decay 
in these minimal SO (10) models, gotten by imposing the coupling constant unification con- 
dition. First we show that, contrary to what is found in the current literature 0, having a 
mixing via the minimal 2x2 mass matrix among the color triplet Higgses does not of itself 
give rise to a new degree of freedom which can be varied in order suppress the proton decay 
rate. Second, we show that the intermediate (~ 10 10 GeV) mass states which often arise 
in these models, coming from weakly coupling the two gauge breaking sectors, leads to an 
enhancement of proton decay rate. These results, combined with the fact that proton decay 
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is already enhanced considerably in SO(10) models due to the largeness of tan (3, pushes the 
values of the parameters of these models into unnatural ranges when the experimental bound 
on proton decay is taken into account. Below, we will demonstrate these points within a 
general framework. We will show elsewhere a more explicit demonstration of these points 
using the specific model by Babu and Barr [Q. 

First we consider proton decay. In minimal SO(10) models, there are two 10's of Higgses, 
each of which decomposes under SU(5) as 5 + 5. Generally, if there is more than one 5 
+ 5 pair, there is a possibility of mixing between the triplet partners of the light doublets, 
that couple to the quarks and leptons, via a mass matrix. When such a mixing occurs, 
the superpotential may be written schematically as Hi J + KHi + AijjHjHj, where the 
5+5 pairs are denoted by the indices (/, J), and we have chosen a basis where Hi and H 1 
are the Higgs which couple to the quarks and leptons. The quark and lepton currents are 
denoted by J and K and M. is the mass matrix. Integrating out Hi then gives — K{M.^ 1 )u J. 
Using a somewhat less bulky notation, 1/Mpo = (.M _1 )ii, we then find two baryon number 
violating dimension four terms in the superpotential suppressed by Mpo, identical to the 
usual minimal SU(5) case with the mass of the color triplet, M# 3 , replaced with Mpj> 

The actual calculation to estimate the proton decay rate is quite complicated and will not 
be repeated here. An extensive calculation was presented in f?]f|, where the p — > i/K + decay 
rate was given by T(p — > ViK + ) = C((3 p / ' M-pT>) 2 \ABi\ 2 . Here, C and [3 P depend upon strong 
interaction parameters, while A depends on electroweak ones. The quantity Bi represents 
the chargino-squark-squark triangle dressing loop, and includes the SUSY mass dependences. 
We note here that there is a sin 2/3 factor in its denominator. Thus, the the decay rate is 
proportional to 1/(Mpd sin(2/3)) 2 which, for large tan/?, is ~ (f|^) 2 - An analysis, running 
over the SUSY mass parameter space gives, from Kamiokande data, Mppj > 1.2 x 10 16 GeV 
where the bound involves a lower bound of tan/? to be ~ 1.4 || which comes from the top 
Yukawa developing a Landau pole at the GUT scale. Noting the tan f3 dependency in the 
decay rate, we then require 
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M PD > tan/3 x 0.57 x 10 lb GeV. 



(1) 



Since one expects M PD to be of order of the GUT scale, ~ 2.4 x 10 16 GeV, this bound 
poses a serious problem to SO(10) models for which tan (3 is expected to be O(50), i.e., one 
must arrange to have Mpd significantly larger than the GUT scale. 

In the models we are studying, the Dimopoulos-Wilczek (DW) mechanism is used to 
create the doublet-triplet splitting. The basic idea is to achieve a breaking of SO(10) via 
growing a VEV for a 45 of the form: (A) = diag(a, a, a, 0, 0) <8> ia 2 where A is the 45 
in the 10 x 10 anti-symmetric tensor notation. This VEV breaks SO(10) down to the 
SU(3)xSU(2)xU(l)xU(l) subgroup where the extra U(l) is identified with the B-L sym- 
metry. The extra U(l) is broken in minimal models via VEV growth for a 16+16 pair. In 
the one-step unification picture, we would expect a to be O(Mgut)- 

If we now couple the 45 to a pair of 10's in a natural way and have a mass term only 
for one of the 10's, then the superpotential will be Wh = M 2 H$ + AHiAH 2 where i7's are 
the 10's. Substituting in the VEV for A, the mass terms for the doublet parts (indicated 
by superscript d) and the triplet parts (indicated by superscript t) are then 



Here the bar in the H indicate that it is part of the 5 under the SU(5) decomposition, 
10=5+5. Thus one of the doublets is automatically massless. These mass matrices also 
suggest the possibility of suppressing proton decay by adjusting the parameters A and M 2 . 
The massive parameter which controls Higgsino-mediated proton decay is given by 




(2) 



and 




(3) 




(4) 



4 



where Ai is the triplet mass matrix. By choosing M2 to be smaller than Aa (expected to be 
O(Mgut)); it appears possible to adjust Mpo so that it is significantly larger than Mqut- 
However since these parameters are also correlated to the mass spectrum of the model, which 
affect the coupling constant running, we must check if this suppression does not ruin the 
coupling unification picture. This is the task we now turn to. 

In order to discuss the running of the coupling constants, it is convenient to introduce 
the variables: yi = l/a* = Ait/g? and x = ^ln(///Mz). Then the RGE's ( = /3j) can 
be written, to two-loop order, as 



where the bi are the one-loop (3 coefficient and bij are the two-loop coefficients. 

The one-loop solutions are just straight lines in the x-y plane with the slopes given by 
the one-loop coefficient. The thresholds corrections in mass-independent renormalization 
schemes are then just a change in the slope at the corresponding x's. The (3 coefficients for 
arbitrary gauge theories may be found in [^]. 

The RGE's can not be solved exactly at two-loop, but the solution can be approximated 
by iterating the one- loop result. The error introduced in this approximation is then expected 
to be the size of the next-order correction. The threshold effects at two-loop are expected to 
be small. In the numerical results below, we have considered it in two situations: (1) when 
the threshold involves many particles so that the change in the two-loop (3 coefficients is 
large, and (2) when the threshold involves a large hierarchy in the energy scales (i.e. when 
there exists an intermediate scale). 

The condition for unification can be stated in the DR scheme as simply the coming 
together of the coupling constants to one point |[LG|| . The threshold corrections would then 
make the running of all the coupling constants the same. Calling the running coupling 
constant in the GUT to be ?/gut, the unification condition is then |/gut(^) = yi(%), x > xjj, 
where x\j is the highest threshold. Adding in the threshold corrections at the one-loop level 
can be done with just simple analytical geometry, and we arrive at 



dx 



y 3 



(5) 
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2/gutO) = Sfc(O) + 6gut3; + Syf - b aiX a , (6) 

a 

where 5yf denotes the two-loop level contribution, yi(0) are the inverses of the coupling 
constants measured at M z , &gut is the one-loop (3 coefficient in the unbroken GUT, and 
x a = ln(^) where M a 's are the threshold masses. 

Eq. |6] represents three constraints on the parameters of a GUT model. Of these, we 
always need to assign one to solve for t/gut whose experimental constraint is less stringent 
than the constraint we can impose on it by grand unification. Thus there are two constraints 
upon the remaining parameters of a model. We can describe these constraints in terms of 
three-dimensional vectors which we may dot into the unification condition. It turns out, 
from the study of (3 coefficients, that the two useful vectors are (1,-3,2) and (—5,3,2), 
where the first one is relevant for Mp^. 

We consider first the minimal SU(5) model, the superheavy particles are the ((2,3)_5/ 3 
+ h.c.) components of SU(5) gauge fields (vector multiplet) in the 24 representation (V), 
the ((3,l)o + (l,8)o) components of the Higgs in the 24 (A), and the color triplet Higgs 
((l,3)_2/3 + h.c.) of the Higgs in the (5+5) (H). We have used the notation (SU(2) rep., 
SU(3) rep.)Hypercharge to show the representations under the Standard Model gauge group. 
The three superheavy threshold masses are labeled My, Ma, and M Hi . The proton decay 
mass parameter, Mp^ is here just M# 3 . Dotting the (1,-3,2) vector into the unification 
condition, Eq. || yields the necessary condition for Mpr>: 

( M Hs \ _ (Mpd\ _ /Msusy\5/6 

[ M%y ^mpJ [ m z ) [n 

where M^ is a quantity which is determined mainly by the low energy measurement of the 
coupling constant, i.e. it is what the mass of the Higgs triplet would be for minimal SU(5) 
when SUSY threshold corrections are neglected (Msusy = Mz)- It can be determined by: 

= 777^(6 sin 2 (M-l) ' 



2tt M z Aa(M z y v ' ' Qa 3 (M z ) 

+^(1.-3,2)^- (8) 



This quantity is plotted in Fig. [I] in the sm 2 (9w) — oi^{Mz) plane, using the estimated two- 
loop correction by setting the GUT scale to be Mpo = M# 3 . (We have used a(Mz) = 
1/127.9.) The right hand side of Eq. (|7|) is not expected to be large, i.e. 0(1), and thus 
Fig. [1] may be viewed as bounds on M# 3 as a function of sm 2 (9w) and a^{Mz)- The proton 
decay constraint of Mu 3 > 1.2 x 10 16 GeV (for the minimal SU(5), i.e. small tan/3) indicated 
by the dotted line can be seen to push the satisfactory a 3 (M^) and sin 2 ^^) values to the 
corner of the experimentally allowed region. This is in agreement with the well-known result 
that it is difficult to achieve low values of a^iMz) in the minimal SU(5) model [JTT], 12|| . 



We are interested in obtaining a similar restriction on Mpd for SO(10). For the models 
with two 10's where the doublet-triplet splitting is done by DW mechanism, we have seen 
that M PD = (A 2 a 2 )/M 2 (Eq. |). From the mass matrices in Eqs. ^| and |3], we see that 
M PD is the product of the triplet masses eigenstates divided by the heavy doublet mass. 
Or, using T 3 and H 3 to denote the color triplet eigenstates, and T 2 to denote the heavy 
doublet eigenstate, Mpo = Mt 3 Mh 3 /Mj< 2 . This turns out to be the same combination of 
the threshold masses which appear in the unification condition gotten by using the vector 
(1, —3, 2). Neglecting for the moment the threshold effects due to the splittings among other 
superheavy particles, we find that the analogue of Eq. [7| is: 



( M H , \ _ / M SUS Y \ 5Ay Mr 2 x 

^M§y { M z ) \M n )' [ ' 

or, 

/ M PD \ _ / Msusy \ 5/6 

l M (o)J-l Mz ) ' UUj 

which is exactly the same as in the minimal SU(5) case. Thus, having a DW derived mass 
matrix mixing does not offer an advantage over the minimal SU(5) case in suppressing 
proton decay rate since it is precisely the same combination of parameters which affects 
proton decay that is constrained by the unification condition. 

A loosening of the constraint on Mpd can come from threshold effects arising from the 
splittings amongst other superheavy particles in a model. In general, minimal SO(10) models 



contain vector multiplets in the adjoint (45) representation, plus chiral multiplets in 45, 54, 
16 + 16, 10, and singlet representations. Upon breaking to the standard model gauge 
group, the massive components that can affect the running of the coupling constants are: 
the vector multiplets in ( (2,3)_ 5/3 + h.c. ) [V], ( (1,1) 2 + h.c. ) [Vx], ( (l,3)_ 4/3 + h.c. 
) [V3], and ( (2,3)1/3 + h.c. ) [Vq] representations; and the chiral multiplets in (3,l)o [A3], 



(l,8)o [As], ( (3,1) 2 + h.c. ) [D 3 ], ( (2,3) 1/3 + h.c. ) [D 6 ], ( (l,6)_ 4/3 + h.c. ) [D' 6 ], ( 
(1,1) 2 + h.c. ) [Si], ( (l,3)_ 4 /3 + h.c. ) [B 3 ], ( (2,3)a/3 + h.c. ) [B 6 ], ( (2,1)! + h.c. ) 



[C2], and ( (l,3)_2/3 + h.c. ) [C 3 ] representations where the labels we will use for them are 
indicated within the square braces. In addition, there are two triplets and a doublet (H 3 , 
T 3 , T 2 ) coming from the two 10's. The A's come from the 45, the .D's from the 54, the B's 
are a combination of the 45 and the 16 + 16, and the C's are from the 16 + 16. When 
the unification condition is imposed, we find then that: 



In the one-step unification picture, the threshold corrections are expected to be small. 
In certain SO(10) models, on the other hand, the necessity to not badly destabilize the DW 
form of VEV is satisfied by weakly coupling the two gauge breaking sectors which in turn 
leads to introduction of of intermediate mass scales In these models, the massive states 
B 3 and B e are of intermediate scale while all the other states (including B{) are at the GUT 
scale. Thus for these models, the threshold effect in fact enhances the proton decay rate. 

Now we turn to the possible solutions to the problems raised here. The results above 
point out how the proton decay rate is intimately related to the threshold effects in these 
SO (10) models. Bringing these models into agreement with experiment hence require large 
threshold effects. One is then faced with having to explain two large effects: (1) why the 
value of Mpd = Mt 3 Mh 3 /Mt 2 is much larger than Mgut, and (2) why the threshold effects 
(the right hand side of Eq. (|TTJ) is also large, in the right direction, and of the correct size 
to explain the ratio M-p^/Mjj . 




(11) 
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One possible way to explain both of these effects is simply to adjust the values of the 
threshold masses |13[]. If we take for granted that item (1) is solved by adjusting Mpo to 
be large enough to satisfy Eq. [L] for tan (3 ~ 50, then item (2) can be accomplished by a 
threshold effect to produce a shift in a 3 (M^) > 13%, or > 5 std. as can be seen from Fig. 1. 
If there are intermediate mass states as discussed above, it provides a threshold effect in the 
wrong direction, and the remnant threshold effect must then be > 40%, or > 16 std. in order 
to bring the ratio Mpd/M^ into agreement with experiment. These are to be contrasted 
with the "successful" prediction of a^{Mz) (within < 4% or < 2 std.) within the framework 
of supersymmetric grand unification where thresholds are neglected. Thus the naturalness 
of grand unification in such models is eroded. 

Perhaps a more natural solution would be to introduce a mechanism which suppresses 
proton decay in a different way, so that it is not strongly dependent on threshold masses. 
For example, with three 10's, three 45's and two 54's, Babu and Barr have shown it is 
possible to strongly suppress proton decay |b|]. Such models then require a much larger 
particle spectrum. 

In conclusion, we have shown that the proton decay constraint plays an important role 
for minimal SO(10) models where the largeness of tan/3 requires that there be a mechanism 
to suppress the Higgsino mediated proton decay. We have shown that the mass matrix 
mixing among the two pairs of superheavy triplets that occur naturally in the Dimopoulos- 
Wilczek mechanism does not in itself give rise to such a suppression, since the proton decay 
rate is constrained by the unification condition in exactly the same manner as without such 
a mixing. We have also shown that in those models where intermediate mass states exist, 
proton decay is in fact further enhanced, making the proton decay constraint even more 
difficult to satisfy. 

This work was supported in part by the National Science Foundation under grant number 
PHY-9411543. 
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FIG. 1. Contour Plot of AfJ§ in sin (9w)—a^(Mz) plane. The dashed line corresponds to the 
lower bound on Mpd in minimal SU(5) of 1.2 x 10 16 GeV. The dot-dashed line corresponds to 
the lower bound on Mpd in minimal SO(10) of 2.7 x 10 17 GeV. Currently, the measurements are 
sm 2 (0 w ) = 0.2313 ± 0.0003 and a 3 (M z ) = 0.119 ± 0.003 pi. 
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